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Abstract 

The characterization and evolution of a Spin Density Wave into the Quantum 

Neel ground state is considered in the context of a weak couphng theory of 

the half-filled Hubbard model. Magnetic properties obtained from this weak 

coupling approach in one dimension compare favorably with exact results from 

Bethe ansatz (B A) . A study of the evolution of several length scales from weak 

to strong coupling is also presented. 
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I. INTRODUCTION 



Interest in low-dimensional Quantum Antiferromagnetism experienced a sharp rise with 
the realization that magnetic fluctuations are at the root of many of the exotic properties of 
High Temperature Superconductors. It is fairly well established by now that the plain two- 
dimensional Quantum Heisenberg model can account for most of the experimental features 
of undoped cupratejiHl. It is also widely believed that the simple one-band Hubbard Model 
(HM) can map large parts of the experimental phase diagramBSH^, and most of the theoretical 
work has used it as an appropriate starting pointHi. Simple as it looks, the HM displays a 
rich variety of regimes, with highly non-trivial physics. In fact, there is no consensus about 
the adequate low-energy effective action for small doping concentrationJzHi, nor about the 
gross features of the phase diagram in that regime000. 

The situation is much more clear exactly at half-filling, because charge degrees of freedom 
are frozen out. In this case and for large Coulomb repulsion, the physical electrons are 
localized, and the ground state is the quantum analog of the Neel state. In the opposite 
limit of small U, the ground state is a Spin Density Wave (SDW), and should be adequately 
described using standard RPA theory over a broken symmetry. Mean Field SDW state. In 
an early development, Schrieffer, Wen and ZhangH proposed that the later, weak coupling, 
approach (to be called RPA-SDW hereafter) could also explain the physics of the Quantum 
Neel state. Later on, Chubukov and FrenkelEl explicitly showed that many of the known 
results for the Heisenberg model could indeed be recovered within the SDW-RPA approach. 
There are, nevertheless, some details missing in the physical picture of the half-filled HM. For 
instance: it is common lore that a SDW is the analog of the Neel state when the electrons are 
itinerant, but a clear and quantitative description of it is lacking. As a direct consequence, 
the crossover regime between both states is not fully understood. The author also feels that 
there is some confusion about the Temperature and length scales generated in the different 
regimes. To be specific: these scales are, in general, well established for the positive-?/ HM 
at strong coupling and for the negative-^/ case at weak coupling, but not in the other two 
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Iimits0. 

The purpose of this paper is to shed some more hght on these issues. In order to do so, 
I will apply the RPA-SDW expansion used in ref. 12 to the one-dimensional case at Zero 
Temperature and compare with the Bethe ansatz solution^ Hence, the whole phase diagram 
at half-filling will be covered. I will also draw many consequences on the scales involved in 
the problem by the heavy use of duality relations^. The computations are performed for 
the worst case, because RPA schemes are supposed to do a bad job in low dimensions and 
in the spin-disordered phase. But I will actually show that RPA works remarkably well even 
in Id -at half filling-. 

Some comments are in order at this point of the introduction. The HM at half-filling 
falls into the 0(3), instead of XY, universality class. A closely related fact is that its 
ground state is a singlet in the strong coupling limit0 so that all correlation functions are 
rotationally invariant. A further, specific feature of one-dimensional systems is that they do 
not possess Long Range Order, even at zero Temperature, because of thermal and quantum 
fiuctuations0. Correlation functions related with the order parameter have a mass gap and 
do not follow Goldstone behavioJil. The complications derived from spin-charge separation 
are eluded because charge degrees of freedom are frozen out at half-fillingEl. 

The layout of this paper is as follows: several magnitudes will be computed within RPA- 
SDW and compared with Bethe ansatz results in section 2. Section 3 will present a study of 
the different length and Temperature scales of the problem. A brief conclusion will end the 
paper. Details of the calculations will be relegated to appendices A and B. Energies will be 
measured in units of the tunneling amplitude, t. 

II. STUDY OF THE SDW-RPA ANSATZ AT T = 

It is commonly believed that RPA approaches should not perform well in low dimensions 
or when one is trying to describe physical properties of the strong coupling phase. The 
purpose of this section is to show that RPA-SDW produces accurate results at either weak 
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or strong coupling even in the extreme case of one dimension, as long as the ground state is 
close enough to a commensurate Antiferromagnet. 

The first quantity which arises in RPA-SDW theory is the magnitude of the Mean Field 
charge gap, 2 A, which is obtained from the equation 

-=-r- (1) 

Inclusion of RPA fluctuations does not modify the equation. Although the asymptotic 
expression for A when U approaches zero, 

sjj^~2n/u^ is not captured by the RPA-SDW 
approximation, which gives A ~ ^-2^/^^ ^j-^g numerical differences are negligible for such 
small values of U (see fig. 1). The RPA-SDW curve begins to depart appreciably from the 
exact result for values of U of the order of the bandwidth, D. This fact suggests that the 
crossover from nearly itinerant to localized electron behavior is not described quantitatively 
within this weak-coupling approximation. 

The exact large-f/ value of the charge gap, on the other hand, is U — At = U — D, 
where D is the free band width@. RPA-SDW theory overestimates 2 A by because the 
approximate lower and upper Hubbard bands tend to (5-fuctions while the exact ones retain 
a bandwidth of D when U ^ oo (see fig. 1 again). 

The collective excitations of the system appear as poles of the imaginary part of the 
different correlation functions x"('?;'M^)- The weight of a particular fluctuation {q,w) on the 
ground state is given by the residue of x" at that particular value of (g, w). The RPA-SDW 
ansatz, in conjunction with Goldstone's theorem, tells us that the only low-energy collective 
excitations of the half-filled HM at zero Temperature and for large U are spin waves, and 
show up as poles of the transverse spin susceptibility, Xrpa" ■ Charge and longitudinal 
spin fluctuations correspond to poles of Xrpa" Xrpa \ have a mass gap. These 
correlation functions have a very simple form in the strong coupling limit. For instance 

1 ^^-%^ + 4(l-cos(g)) [ 1 

Fig. 2 (a) shows the imaginary part of the transverse spin susceptibility for U = 20. It 
displays the typical Goldstone-mode structure and has a strong q-dependence, proper for 
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spin waves. Notice also that the weight of the Ferromagnetic zero mode tends to zero, while 
that of the Antiferromagnetic one diverges, in agreement with strong coupling analyses^. 
On the contrary, Xrpa" ^ = 20 is almost dispersionless and has a large gap of 

order U (see figure 2 (b)). Hence, longitudinal spin fluctuations are massive and possess a 
localized nature. 

The exact spin correlation functions should be rotationally invariant because the ground 
state of the system is a singlet at strong coupling. Regretfully, the RPA-SDW response 
functions do not have this property. 

Returning to the strong-coupling case again, it is easy to show that the excitation spec- 
trum given by RPA-SDW is 



This is the result one would obtain for the Heisenberg model in Linear Spin Wave theory 



and serves to define Vg for all values of the Coulomb repulsion (fig. 3). Notice that the 
RPA-SDW velocity grows as U increases from zero and then has a maximum for [/ ~ 1. 
What happens here is that the limit U ^ 0, where Vg = 2, does not coincide with the U = 
case {vs = 0). This subtle point is well captured by the exact BA analysis, but not by 
RPA-SDW theory, where Vg begins at zero, increases steeply until it reaches a value close to 
3 at the maximum and then approaches the BA curve for for U ^ D. Both curves decrease 
approximately at the same pace for even larger U and tend asymptotically to |J ~ 1.6 J 
and J. The value given by RPA-SDW theory for Vs in higher dimensions {\plJ for d=2) 
coincides with that obtained from LSWT. 

The squared magnetic moment, < S*^ > is displayed in fig. 4, as given by Mean Field 
Theory (dashed line) and Bethe ansatz (solid line). It is quite remarkable the fact that Mean 
Field Theory already give the right f/ — > and f/ — > oo limits. Notice also how Mean Field 
Theory underestimates grossly < S*^ > for small values of [/, while tends correctly to 3/4 
in the opposite limit. What this seems to imply again is that the Mean Field spin degrees 



Wq = Vs\ sin(g)| ~,^o v^q 



(3) 




I. The linear relation between w and q for small q holds all the way down to f/ = 
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of freedom are too itinerant in the crossover regime from itinerant to locahzed electron 
behavior. 

RPA fluctuations do not cure this flaw. In fact, < S"^ > diverges at strong coupling as 

o2 1 / 1 - cos(fc) 
< >= -} — . oo (4) 

which is also the result obtained by LSWT. In the opposite limit, on the contrary, RPA-SDW 
theory reproduces the exact free result |. 

Another important result concerns the average of the staggered magnetization for large 
U, which also coincides with LSWT0 

<S'>=-il-2j:''-f^]^-oo (5) 
2\ V \Mk)\ J ^ ^ 

This result implies that one can obtain a charge insulating spin disordered ground state, 
also called spin liquid, within SDW-RPA theory, because the charge gap is still given by the 
Mean Field result. LSWT is known to diverge in one dimension due to infrared divergences, 
but to lead to amazingly accurate results in higher dimensions@'0. 

Figure 5 is a useful eye-guide to understand the differences between the Neel state and a 
Spin Density Wave. For large values of U, the charge degrees of freedom of physical electrons 
are frozen out and the spin degrees of freedom become localized magnetic moments: the 
Heisenberg Hamiltonian provides a good description of the HM in this limit. In this case, 
< >= S{S + 1) = 3/4. In the opposite case, < S*^ > equals 3/8, which is the value for 
free electrons. As soon as U differs from zero, a charge gap opens up, and the expectation 
value of the staggered magnetization -in dimensions higher than one- becomes finite albeit 
small. The squared magnetic moment increases with U from 3/8 to 3/4 continuously but 
with a steep slope. By the time when U equals the band width, < S*^ > is already very 
close to 3/4, and the Spin Density Wave has fully developed into the Neel state. Likewise, 
for that value of U, the exact spin wave velocity almost coincides with | J. 
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III. LENGTH SCALES IN THE HUBBARD MODEL 

The transformation 

/it ^ fi\ 

ft ^ e^^^7+ (6) 

where Q = (vr, tt, . . .), changes the positive-U half filled HM into the negative one, and vice 
versa. It also exchanges spin and isospin operatorJl^: 

<i,j>,a i <i,j>,cr i 



These duality relations allow to infer many properties of one model from the other. For 
instance, the positive-?/ HM at half-filling possesses a global SU{2)sp ® U{l)ch- This means 
that its counterpart has a global SU{2)ch ® U{l)sp- Likewise, an 

• Antiferromagnet (AFM) polarized along the Z-axis is equivalent to a commensurate 
Charge Density Wave 

• AFM polarized along the XY-plane is equivalent to a uniform superconductor 

• Ferromagnet (FM) polarized along the Z-axis is equivalent to the uniform state 

• FM polarized along the XY-plane is equivalent to a commensurate staggered 
superconductor^ 
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Because the ground state of an AFM at strong coupling is a singlet, a uniform super- 
conductor should be degenerate with a CDW for large U. Conversely, the spin transport 
properties of an AFM polarized along the Z-axis and XY-plane should be very different, 
because the charge transport properties of a uniform superconductor and a commensurate 
CDW indeed are. The appearance of the nesting wavevector Q in the duality relations 
introduces, however, a slight difference because staggering means folding of the Brillouin 
zone into two pieces, with the corresponding doubling of the periodicity. It is also interest- 
ing to notice that there do not seem to be many experimental realizations of an staggered 
superconductor, while a uniform FM polarized in the XY-plane takes place in Nature very 
frequent lyii. 

The remaining of this section is devoted to elucidate what are the different length scales 
which appear in the positive and negative-U HM and how they show up in the following 
correlation functions: 

\ 



GsDw{k, k + Q) 



< TUk)f+ik) > < Tf^{k)f^{k + Q)> 
< Tf^ik + Q)f+{k) > < TUk + Q)t{k + Q)> 



G sup{k^ 



^ Go f] ( < TUk)t{k) > 



F+ Gn 



< Ti){k) > 
< Tij+ik) > < Tt^{k)f+^{k) > 



C+-{k) =<TS^{k)S-{-k) >^C^{k) =<T^+{k)ij{-k)> 



(8) 



The language proper for superconductivity will be used, and it is assumed that the 
dimension is larger than 2. RPA-SDW makes the ansatz that the spin is polarized along 
the Z-axis. However, I have chosen in this section to write down Green-functions which 
correspond to polarization in the XY plane, Eq. 8, to make duality arguments more explicit. 

Generically, there are three different Temperature scales present in this problem, as 
shown schematically in figure 5. The highest Temperature, denoted as T^, is that at which 
binding of uncorrelated pairs occurs. It is a crossover Temperature, signaled by a divergence 
in the T-matrixBil. A spin gap, mi = 2A, opens up in Go because it costs some energy 



8 



to break a pair. Conversely, for U positive, a charge gap appears in Gq. The non-diagonal 
elements of Ggup and Gsdw: F{k) =< T-ip^{k) >, are still zero. 

At a lower temperature, Tq (~ J for strong coupling), correlations among Cooper pairs 
build up, and the correlation length ^(T) begins to be longer than the lattice spacing, a. 
Its inverse, 1171^0^.^ = 1/^, appears as a pole in the correlation function C^, which takes the 
conventional Ornstein-Zernike form 

(9) 

nicorr is the gap for excitations of the collective modes of the order parameter. For U positive, 
the pole appears in the correlation function C"*"~, and nicorr is the gap in the spectrum of 
spin wavesS. This is also the Temperature where the anomalous Nambu Green's functions, 
F{k ~ 1/a) =< Til)^{k ~ 1/a) > begin to differ from zero. 

Finally, at a still lower temperature, Tc, a phase transition towards the broken symmetry 
superconduting state occurs, and is signaled by a divergence of ^(T): the mass gap nicorr 
closes down, as dictated by Goldstone's theorem. Below Tc, another length scale sets in: the 
coherence length of Cooper pairs, which is also Josephson's coherence length, = 1/mji^: 

C.,„).4.0^.?^^ (10) 

where if)* =< il)^{k = 0) > is the expectation value of the order parameter and ps is the 
superfiuid density. 

Notice that for small U , T^, Tq and Tc merge into one line; rrii is then called the supercon- 
ducting charge gap, and equals m j for T < Tc, or nicorr for T > Tc. This means that ^(T) 
decreases to zero very quickly as the Temperature is raised above Tc in the weak coupling 
regime. In other words. Cooper pairs become completely uncorrelated for Temperatures 
slightly above T^. Mean Field BCS theory is very accurate in this limit and predicts that 
the order parameter = 

For larger attraction U, these three Temperatures depart from each other, and there is a 
finite Temperature range above Tc where ^(T) is different from zero. The order parameter 
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is no longer equal to ^ because fluctuations become important (remember the discussion 
in the former section). 

Finally, for very large U, diverges, and Cooper pairs become composite bosons. Tq 
and Tc again go to zero because it costs more and more energy for a pair to hop from one 
site to the following -the effective tunneling amplitude, t^/f/, vanishes-. For the positive-U 
case, we would say that the spins become progressively less correlated because the exchange 
integral, J ~ 4t'^/U vanishes. 



IV. CONCLUSIONS 

This paper shows that the RPA-SDW expansion developed in ref. 12 provides reliable 
results for the magnetic properties of the Hubbard model at either weak or strong coupling 
even in the extreme case of one spatial dimension. RPA-SDW does provide a qualitative, 
yet not quantitative, description of the intermediate coupling regime, when the Coulomb 
repulsion U is of the order of the bandwidth D. These conclusions should be quite robust 
against doping and Temperature as long as the exact ground state is still Neel-like000. 
The accuracy of this expansion should improve with increasing dimensionality. The Spin 
Density Wave and Quantum Neel ground states are characterized the expectation value 
< S"^ >, which measures the degree of itinerancy of the charge degree of freedom electrons. 
This criterion allows to show that the Spin Density Wave evolves continually into the Neel 
state. Free electrons have < S*^ >= 3/8, but < S*^ > increases quickly with U and saturates 
to the value of localized spins, 3/4, when the Coulomb energy U is of the order of the 
bandwidth. Several Temperature and length scales are discussed in both weak and strong 
coupling limits. 

It is a pleasure to thank G. Gomez-Santos, N. Lorente, J. Perez-Conde, F. Gebhard, F. 
Sols and E. Miranda for their helpful hints. This work was supported by DGICYT, Project 
No. PB93-1248 
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APPENDIX A: REMINDER OF SDW-RPA FORMALISM AT HALF-FILLING 

AND ZERO TEMPERATURE 



The Hubbard Hamiltonian can be rewritten in terms of the staggered magnetization, rhi 

n^Yl ^k,M,aka -uY^ml ^n^-^tY. cos(A;)/+ J,,. (ai) 

k i k,cr i 

if we neglect charge fluctuations, a good assumption for the half-filled case. The staggered 
magnetization can, in turn, be expressed in terms of spin- 1/2 operators 

— (Pi) ^t-i ^i ) = X/ fhOi^apfh/B 

= (s.p, ^) (A2) 

as rhi = e'^^^'iSf. A Mean Field decomposition on the interaction term is performed now, 
771^ —< rhi > TTT'i, so that 

^ = E ^kfk,Jk,a + U <rhi>Y^rhi (A3) 

k,a i 

where U < rrij >= A is half the charge gap. The Hamiltonian is diagonalized with a 
conventional Bogoliubov transformation and the result is 

k k,a 

E, = V4 + (A4) 

with the self- consistency condition 

The prime in the sum means summation in the restricted magnetic Brillouin zone. The 
new ground state, \SVyV >, is annihilated by c and ci"*". The self-consistent equation can 
be solved exphcitly for large U and gives A = C//2, which implies that the Mean Field 
staggered magnetization is 1/2 in that limit. 
We define Matsubara correlation functions as 
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X'''{q.<i.r) = -<TS''{q,T) S-'{-q',0) > 



(A6) 



It is straightforward to obtain from them the retarded non-interacting response functions 

1 1 1 



k V ^ki^k+q J 



w + Ek+q + Ek + iS w - Ek+q - Ek + iS ^ 



^ k \ ^k-t^k+q J 



1 



1 



w + Ek+q + Ek + iS w - Ek+q - Ek + iS ^ 



Xq {q,q',w)^ 



Sq,q'+Q ^Y^ffj__^ 1 



-7 + 



2 ^ \Ek ' Ek+q J {w + Ek+q + Ek + iS w - Ek+q - Ek + i5 _ 



(A7) 



The RPA-SDW expressions are now given by 



X^RPA 



1 + f 



^,zz Xo 

Xrpa 



i-2c/xr 



Xrpa 



Xo (Q:'^) Xq (q^w) 
^Xq~{q,w) Xo~{<1 + Q,w) j 



\ 



(A8) 



-Uxq (q^w) i-Uxo {Q + Q:w)j 

The poles of Xrpa give the spin- wave excitation spectrum in the RPA-SDW approxima- 
tion. They are located around the antiferromagnetic, go — and ferromagnetic, go = 0, 
wave- vectors. Both have the same velocity, but different weights. Expanding the denomina- 
tor for small g' = go — tt and the spin- wave velocity, Vs, in the g' — > limit can be written 
as 



= 4 



lyzU 



X 



(A9) 



where x, y and z are given by 
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1 

k 



^ El 



^,3sin^(/c) - 1 ^,cos^(A;) sin^ffc) 

y^2^ — 12 }^ -5 

k ^k k ^k 

^ = E'^ (AlO) 



This expression for the spin-wave velocity gives the correct asymptotics for a half-filled band, 

— VdJ U ^ oo 
t 

^0 U ^0 

where J = ^ is the exchange integral. 

The summations over the restricted Brillouin zone can be performed analytically for 
large U. A straightforward yet tedious calculation gives the following expression for the 
transverse correlation function 

+ - f \ 1 ^/;^-< + 4(l-cos(g)) [ 1 1 \ 

U 2Wq [W + Wq + ld W — Wq + id ) 

where Wq = Vs \ sin(g)| is the spectrum of spin waves. The imaginary part of the transverse 
response function gives the spectral density of these low-energy excitations: 

X W (Q, tan(|)5(«; - Wq) (A12) 

Notice that Xrpa' (?) '"^) ^ ^^^^ mode at q — with vanishing weight and another at 
q — TT, whose weight diverges. The real part of the response function at zero energy 

XW(?,0)^^1 (A13) 

is related to the transverse correlation function C+ ~ . 

The RPA-SDW staggered magnetization can be written as 

<S'>^<mi>^^J2a [ —A^{k,k + Q,w)nFiw) = ^J2[ ^A^{k,k + Q,w) 

(A14) 
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where the spectral function Acr{k, k', w) is obtained from Dyson's equation. The RPA-SDW 
self-energy is given by the matrix equation 

t^{k, w) = ^Y. Gl{k + k\w + wb) X-.{k', wb) (A15) 

^ k',Wg 

where only the transverse and longitudinal spin plus charge fluctuation channels are taken 
into accounts. It is important to realize that RPA fluctuations reduce the staggered mag- 
netization from its Mean Field value, < >= A/U. In particular, for large values of U, 
Mean Field Theory yields < 5^ >= 1/2 while RPA-SDW gives 

This result implies that one can have a charge insulating spin disordered state, or spin liquid, 
within SDW-RPA theory, i.e.: < S"^ >< even though A ^ 0. 

Finally, the expectation value of < S*^ > is defined in terms of the dynamic correlation 
function as 

<s'>-<s^>'=j:r ^ c{q, w) (A17) 

which, in turn, is related with the dynamic response function through the fluctuation- 
dissipation theorem which reads at T = 0: 

<S^>-<S^>=Y. —{X'''''{q,w) + x'-''{q,w)) (A18) 

q JO 

The double prime denotes again the imaginary part of the function and the summation is 
now performed over the whole Brillouin zone. It is interesting to notice that the Mean Field 
approach gives both the weak- and strong-couphng limits correctly 

= + ^ + ^ U = oo 
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APPENDIX B: BETHE ANSATZ RESULTS 

The Hubbard model was solved by Lieb and Wu (LW)§ in 1968. Since then, a number of 
authors have completed LW's work and obtained the different physical magnitudesH'&il. 
This appendix is devoted to display the integral equations, the charge energy gap, < S*^ > 
and Vg as obtained by Shiba0 and OvchinnikovS. 

The integral equations at zero field and half-filling for the distribution of the k-rapidities 
are given by 

p{k) = — + - cos{k) r dk' p{k') S{2{sm{k) - sin(A;'), u) (Bl) 
duPik) = — cos{k) / dk' p{k') _R(2(sin(fc) — sin(A;'),-u) + 

TIU J-n 

+ - cos{k) r dk'dup{k') {S{2{sm{k) - sin{k')),u) - S{2siTi(k),u)} 

IT J-n 



where the kernels S and R are 

oo 

5(x,«) = E(-ir 



2 f^„,\2\ 



i.(...)^B-ir' t."':!:!i;r 



[x'' + (nu) 

The charge energy gap is given by 



- = l6u dx ^ (B2) 
t Ji sinh(^) ^ ^ 



while the velocity of spin waves in the q limit has the following expression 

where J„ are Bessel functions of imaginary argument. Finally, 

3 



(B3) 



< >= - + 3 r dk{l + cos(A;)) dup{k) (B4) 

4 J— TT 
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FIGURES 

FIG. 1. Charge gap, 2 A as a function of U. Solid line is the bethe ansatz exact answer; dashed 
line is the RPA-SDW Mean Field result 



FIG. 2. (a) Xrpa ^ ^ function of q and w iov U = 20; (b) X^rpa il^ same value of U 

FIG. 3. Spin-wave velocity, Vg, versus U. Solid line is the exact result; dashed line is RPA's 
answer. Notice that Vg tends to as U goes to 0, and shows a maximum at C/ ~ 1 

FIG. 4. < 5^ > versus U. Solid line is the Bethe ansatz result and dashed line is Mean Field's 
answer 

FIG. 5. Characteristic Temperature of the different phenomena as a function of U. Solid line 
is the Binding Temperature, T^; dashed line is the Temperature, Tcorr, at which correlations begin 
to build up; dotted line is the transition Temperature to Long Range Order, Tc 
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